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Abstract 
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1 Introduction 



Although all experimentally observed particles are either bosons or fermions, the general 
principles of quantum mechanics do not prohibit existence of particles obeying other types 
of statistics 0]. The motivation to consider an exotic statistics is that it may provide an 
effective description of the dynamics of particles if their interaction is so strong that the 
available occupancy for a given state depends on the number of particles already present 
in the state [^. An exotic statistics can also arise if particles possess a hidden internal 
degree of freedom that is invisible in the Hamiltonian but that can become dynamically 
relevant. Finally, an exotic statistics can emerge in description of physical models in one or 
two spatial dimensions. The reason is that in this case the multi-particle wave function is not 
necessarily even or odd and particles (anyons) may obey a fractional statistics Q. Moreover, 
in low-dimensional theories the exchange statistics of the fields present in a Hamiltonian is 
not directly related to the exclusion statistics of the corresponding particles. For instance, 
although particles in the one-dimensional real coupling affine Toda field theories (ATFT) are 
bosons, one has to impose the Pauli principle on their momenta in the thermodynamic Bethe 
ansatz analysis in order to obtain correct values of corresponding central charges |^ . Another 
example is the conformal field theory, where quasi-particle spectra can be constructed with 
use of an exotic exclusion statistics Q. 

The first attempt to introduce an exotic statistics is attributed to G. Gentile Q, who 
proposed an intermediate statistics, in which at most G indistinguishable particles can occupy 
a given state. This intermediate statistics interpolates between fermions and bosons, that 
are recovered for G = 1 and G = oo. A decade after Gentile's work H.S. Green introduced 
the parafermi statistics Q that also fixes maximal occupancy for a given state. Since then, 
various aspects of systems governed by the Gentile or parafermi statistics were discussed in 
the literature, among them possible applications to the elementary particles theory |Q and 
to the quark confinement problem in the QCD 

Many authors discussed thermodynamic properties of an ideal gas obeying the Gentile or 
parafermi statistics |6|, [l^, 11 1. In the present paper we will study the thermodynamic limit 
of one-dimensional relativistic integrable field models governed by a Gentile-like statistics. 
More precisely, the main purpose of this manuscript is to implement systematically a gen- 
eralized extensive exclusion statistics into the thermodynamic Bethe ansatz (TBA) analysis 
of systems where the interaction of particles is relativistic, short-range and characterized by 
a factorizable scattering matrix. 

The TBA was originally developed in the seminal papers by Yang and Yang [12| for 
treating a one-dimensional ideal gas. In later works []l3| |l^ ] this technique was adopted 
to one-dimensional relativistic models. The boundary condition for a many-particle wave 
function leads to what is commonly referred to as the Bethe ansatz equation, that provides 
the quantization condition for possible momenta of the system. The TBA constitutes an 
interface between massive integrable models and conformal field theories. From the TBA 
one can extract information about the ultraviolet limit of the underlying massive model, in 
particular, find the corresponding (effective) central charge. 

In the derivation of the TBA equation the underlying exclusion statistics is usually taken 
to be either of bosonic or fermionic type |^, 12, |l3|. In the present paper the TBA equation 
will be derived and studied for significantly more general situation. Namely, let W{N,n) 
denote the total dimension of the Hilbert space for a system of n indistinguishable particles 
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that can occupy different states in tlie Fock-space. Assume that there exists a function 
such that its iV-th power is a generating function for the dimensions W{N,n), i.e. 

{f{t)f = J2W{N,k)t'. (1) 

A;>0 

In the statistical mechanics, if the variable t is understood as fugacity, the sum on the right 
hand side is the grand partition function of the system. The property of a grand partition 
function to be an A^-th power of a function independent of N is called extensivity. As we are 
going to demonstrate below, this property allows us to develop the TBA analysis for a rather 
general choice of f{t). Let us stress that for the purpose of deriving the TBA equation even 
asymptotical extensivity suffices, i.e., equation should hold in the large N limit. 

Of course, if the explicit form of f{t) is known, we can use it to simplify the general 
formulae. For instance, an instructive example is the Gentile statistics of order G which is 
defined by the following choice of f{t): 

fG{t) = l + t + t^ + ... + t^ . (2) 

For G = 1 and G = oo this statistics describes ordinary fermions and bosons, respectively. 
Although fait) is the simplest non-trivial choice of f{t), it has all features of a generalized 
extensive statistics. Therefore, the Gentile statistics will be our main working example below. 

The paper is organized as follows. In Section 2 we describe possible types of a general- 
ized extensive statistics and study properties of the corresponding entropy densities in the 
thermodynamic limit. The Gentile statistics and the 7-ons statistics are discussed in this 
context. In Section 3 we compare properties of the Gentile statistics and the Haldane-Wu 
statistics. In Section 4 we derive the thermodynamic Bethe ansatz equation and the finite- 
size scaling function for a relativistic multi-particle system in which statistical interaction is 
governed by a generalized extensive statistics. In Section 5 we derive the Y-system related 
to such a generalized statistics. In Section 6 we study the ultraviolet limit of the generalized 
TBA equation and obtain an expression for the corresponding effective central charge. In 
Section 7 we compute finite-size scaling functions and central charges related to some affine 
Toda models for the Gentile statistics and the 7-ons statistics and compare them with their 
counterparts for the Haldane-Wu statistics. Our conclusions are stated in Section 8. 

2 Types of extensive statistics and entropy density 

Consider a system of n indistinguishable particles that can occupy different states in the 
Fock-space. Assume that the particles obey the Gentile statistics of order G, i.e., that each 
state can be occupied by at most G particles (with G being a positive integer number). 
There are several combinatorial ways to compute the total dimension of the Hilbert space 
for such a system. For instance, we first choose mi states which are occupied by at least 
one particle, then among these rrii states we choose m2 < mi states which are occupied by 
at least two particles, etc. This way of counting yields the following expression for the total 
dimension of the Hilbert space 

Ar>mi>...>mG_i>0 
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where = k\/{m\{k — m)l) are binomial coefficients (and CJ^ = if m > /c or if m < 0). 
Now, let us compute A^-th power of the polynomial fait) defined in (^). Thanks to the 
obvious recursive relation fait) = l + t/G_i(t), we can do this by G consecutive applications 
of the binomial formula. Then it is easy to see that A^-th power of fait) is a generating 
function for the dimensions WG{N,n), i.e. 

GN 

ifait)f = J2WGiN,n)e. (4) 

n=0 

Thus, the Gentile statistics is a particular case of an extensive statistics. 

Let us stress that the order G of the Gentile statistics must be a positive integer. Although 
we could try to extend it to other values by formal replacement of (|2|) with 

1 - 

fa{t) = -^^, (5) 

the resulting function fait) will not satisfy the important positivity property (see below). 

Consider now a more general system of n indistinguishable particles that possesses, at 
least for large A^, the extensivity property (||) with 

m = E Pkt' , (6) 

fc=0 

where d can be a finite positive number or infinity. The Taylor coefficients Pk are fractional 
dimensions of levels in a single-state Hilbert subspace and can be regarded as probabilities 
of occupation of the given state by k particles (see related discussions in [|^, [l^, 17]). The 



total dimension of the single-state subspace is /(I). It is natural to require positivity of the 
quantities Pk'- 

Po = l and Pfc > for k>l. (7) 

The first condition here implies that the vacuum is realized with probability one indepen- 
dently of the size A^ of the system. Furthermore, we assume that the series (^ converges in 
the complex plane for < |z| < R, where R < oo. This implies that f{t) belongs to one of 
the following types: 

I. d < oo , R = oo ; 

II. d = oo , R = oo ; , , 

Ilia. d = oo , 1 < R < oo; ^ ^ 

Illb. d = oo, 0< R<1. 

The types I, II, and III consist of finite degree polynomials, analytic functions with infinite 
Taylor series, and meromorphic functions, respectively. For example, a finite order Gentile 
statistics belongs to the type I, the Boltzmann statistics (/(t) = expt) belongs to the type 
II, and the bosonic statistics (/(t) = (1— 1)~^) is of the type Illb. The reason we divided the 
type III into two subtypes will become clear below. 

Let us remark that, in addition to (^), one may require that Pi = l with the motivation 
that the statistics should not be deformed if only one particle is present in the system. 
Although all the concrete cases which we consider below do fulfill this requirement, it does 
not seem to be crucial for the general considerations. 

The entropy of the system under consideration is S = klnW{N,n), where k is the 
Boltzmann constant. We need to compute the entropy in the thermodynamic limit, i.e.. 
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when N ^ oo and n/N is fixed. The extensivity property (||) ahows us to express W{N,n) 
as the fohowing contour integral 

WiN,n) = ^ <f dz {f{z)f = ^ i - exp{iV(ln f{z) - /i In z)} , (9) 

■^■^Ti J\z\=p ^TTZ J\z\=p Z 

where p < R, and we denoted fi = n/N. Now we can find an asymptotic expression for 
the entropy in the thermodynamic limit by application of the saddle point method to the 
integral (^). The saddle point of the exponent in (^), denote it x, is found as the positive 
root of the equation (prime stands for a derivative) 

xf{t)\t=. = ^lf{x). (10) 

Choosing p = x in (0), so that the saddle point belongs to the integration contour, we apply 



the standard result (see, e.g., ||18|) for the Laplace integral: 



W-(iV,MiV) = ^"P^^(;"^(^^-^^""^^ C- + Oil/N)) as N^oo, (11) 

^2TTNh^{x) 



^x 



where hf,{t) = f"{t)/{nf{t)) + (1 - /i)/*^. For any f{t) satisfying (0) and x obeying (|^) 
we have h^{x) >0 due to the Hadamard theorem |l^, |l9|. From (11) we obtain the entropy 
density s{p) for a given value of /i 

lim —\nW{N,pN)=s{ii) = \nf{x)-fi\nx, (12) 



where x is positive and satisfies (|1( 

Two remarks are in order here. Strictly speaking, if f{t) is degenerated, i.e., if there 
exists integer p > 1 such that f{t) = J2k Ppkt^'' ^ then the above derivation modifies since 
f{t) has p saddle points. However, the corresponding physical interpretation is simply that 
particles can be added only in p-tuples. Therefore, it is natural to introduce cluster variables 
t' = pt, n' = pn, etc. In terms of these new variables formulae (|9|)-(|ri|) remain valid. 

More important remark is that for deriving (|l^) the property Pk > is crucial. In general, 
without the positivity property does not have a real and uniform in p limit. Moreover, 
if we allow f{t) to have negative Taylor coefficients, then some W{N,n) in the expansion of 
(/(i))^ may also be negative. This leads to an immediate problem with the definition of the 
entropy as a logarithm of W{N, n). However, if the first negative coefficient in the expansion 
of {f{t))^ appears only at sufficiently large power no (of order 0{N)), then f{t) violating 
the positivity condition can still have well defined entropy density for certain range of fi. 
For instance, such is the case of the Haldane-Wu statistics with statistical interaction g (see 
Section 3), where no ~ N/g. However, in the present paper we will restrict our consideration, 
for simplicity, only to partition function strictly obeying the positivity condition. 

The entropy density is a key quantity for developing the TBA analysis. From the def- 
inition of fj, it is clear that < fi < d for f{t) of the type I and < ^ < cxd for the other 
types. Formally this follows from equation ( [lO| ) which shows that /i varies from /i(x = 0) to 
fj.{x = R). Moreover, /i(x) is a strictly increasing function because dxH = x///i^(x) > for 
x>0. 

Since fit) is regular at zero, it is easy to see that jilnx tends to zero for small fi. Together 
with the condition Po = 1 this yields s(0) = 0. In order to analyze further behaviour of s(/i). 
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we employ (|lO| ) and derive from (|l^ ) that 

d^s{i^) = -Inx. (13) 

Since dxfJ'>0, we conclude that d'^s{fi) <0, that is s(/i) is a convex up function. Furthermore, 
( p!^ ) implies that maxs(;u) = In /(min{l, i?}). Note, however, that s(/x) is not guaranteed to 
be positive for the whole range of fi. More precisely, s{fi) may have one root at the interval 
fi{l) < fi < fJ-iR)- Indeed, for f{t) of the type I we derive from (|lO|) and ( [T^ ) that 

s{d) = In Pd. (14) 

That is, s{fi) is non- negative everywhere at 0<fi<d only if the last Taylor coefficient Pd is 
greater or equal to one. For the other types of f{t) we have the following properties 

s{oo) = oo for type Illb, (15) 
s(oo) = — oo for types II and Ilia. (16) 



Indeed, in the first case x<R<l, hence, according to (|13D, s(/x) is strictly increasing. In the 
second case we have s(oo) = In /(I) — J^^^ dfi Inx. The integral here apparently diverges 
since it exceeds /^ijp)d/i IuRq, where we can choose any Rq such that 1<Rq<R. 

Investigation of thermodynamic properties of models in which the entropy density be- 
comes negative is rather problematic. Therefore eqs. (|l4|) -(|l^) suggest that we should restrict 
the choice of f{t) to the type I with Prf > 1 (the Gentile- like statistics) and to the type Illb 
(the bose-like statistics). However, let us notice that in the TEA framework the variable 
X depends on the rapidity 9. For some relativistic models with factorizable scattering the 
typical picture Q is that x{6) falls off as \9\ grows (for instance, x{9) = exp{— mr cosh^} for 
the trivial S-matrix). In this case x{9) < x{0) and we can consider f{t) of the types II and 
Ilia if s{n) remains non-negative for fi < fi{x{0)). 

In the context of this discussion it is instructive to consider the so-called 7-ons |^, pT| . 
These are particles with statistics interpolating between fermions (7 = 1) and bosons (7 = — 1) 
in the following simple way 

m(x) = . (17) 

I + ^x 

Solving equation ( |To| ) for this choice of fj, it is easy to find the corresponding single-state 
partition function: 

f,{t) = (1 + Jt)'/^ = l + t + ^[^-117^ ' (18) 

k>2 

where [^]7 = nm=i(l~"^7)- -^o'^ 7>0 the positivity condition (|^) is fulfilled only if 7 = l/(i, 
with d positive integer. In this case /'y(t) is of the type I, the corresponding maximal 
occupancy is fi = d. Since InP^ < for d > 1, the entropy density becomes negative for 
certain range of /i. The case of 7 = describes the Boltzmann statistics. Here f-yit) = expt 
belongs to the type II and s(/i) is negative for ^>e. Finally, for 7 < the positivity condition 
^ is always fulfilled and /^(t) is of the type Ilia or Illb depending on the value of 7. 

Returning to the Gentile statistics, we can summarize the properties of the corresponding 
entropy density ^^(/i) as follows (see Fig. 1 for illustration), sdn) is a convex up function 
defined for 0<n<G and vanishing at the end points of this interval. It attains the maximum 
at /i = G/2: 

Sg(G/2) =ln(G + l). (19) 
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Furthermore, Sg(m) possesses the following symmetry 

Sa{G - fi) = sa{^^) . (20) 

This is a consequence of the equality Wg{N,GN — n) = VFg(A^, n) which, in turn, follows 
from the relation (^). 

3 Comparison with Haldane-Wu statistics 

It is interesting to compare thermodynamic properties of a system of particles obeying the 
Gentile statistics with those of a system obeying the so-called Haldane-Wu statistics. For 



the latter, the total dimension of the Hilbert space is postulated to be |21 



s"- ' ^ n]{N - gn + g-l)\ ^ ^ 

This expression interpolates between the bosonic {g = 0) and the fermionic {g = 1) statistics. 
Introduction of such an interpolating statistics was motivated by Haldane's generalization 
Q of the Pauli exclusion principle to the form 

AD/An = -g. (22) 

Here D is the number of available states (holes) before the n-th particle has been added 
to the system. The parameter g is called the statistical interaction. Properties of an ideal 



gas obeying the Haldane-Wu statistics were actively discussed in the literature |15, 21, 22 1 



The corresponding TBA equation for relativistic integrable models was obtained in [23| and 
investigated in |2^, 24 1. 

A direct comparison of the quantities ^^(A^, n) and Wg{N, n) is problematic if the system 
contains a finite number of particles. Indeed, formula (P) has a bulky form for generic G 
(i.e., for 2 < G < oo). Moreover, it requires additional conventions to make exact sense of 
expression (^) for generic g (i.e., for 0<g'< 1). Actually, there exists no prescription for 
counting of states on the microscopical level that would lead to (|2^) (see related discussions 
in ||l^, |l^, 17]). To overcome these difficulties we can compare the two statistics in the 
thermodynamic limit. More precisely, let us compare the entropy density SaifJ-) with its 
counterpart Sg(/i) = limAr_^oo In ( A, /i A) that is readily derived from ( |2T| ) with the 
help of the Stirling formula 

Sgifi) = (1 + - g)) ln(l + ;,(! - 5)) - ^ In /i - (1 - 5/i) ln(l - g^) . (23) 

Let us notice that the Haldane-Wu statistics is extensive but the corresponding function 
fg{t) does not satisfy the positivity property (0). Indeed, we can reconstruct fg{t) first as a 
function of fi by the following formula 

fin) = exp{s(/i) - fi d^s{n)} , (24) 

that follows from (|l^) and (|13|). For Sg(^) given by ( p3D this yields 
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which is a well-defined function. But reexpressing it in terms of t with the help of (|lO|) 
we will always obtain a function that breaks the positivity property (^). For instance, 



for g = 1/2 we find (which coincides with the result of IH) fg{t) = l+t^ /2+t^/T+W/l 
l+t+t^/2+t^/8— 1^/128+. . .. However, the first negative term in expansion of the A^- 
th power of this series appears for n = 2A^ + 3. This implies that in the large limit 
thermodynamic quantities are well-defined for /i < 2. 

For our purposes we can regard (^) as a definition of the Haldane-Wu statistics. Then 
we have a finite maximal occupancy N/ g for a given number of particles N simply because 
Sg{ii) is defined for < ^ < 1/gf. Furthermore, Sg{ii) is a convex up function vanishing at 
the end points of this interval. Thus, we see that properties of the entropy density 
are similar to those of a type I extensive statistics and, specifically, to those of the Gentile 
statistics. Therefore, it is natural to ask how much the behaviour of sdp) differs from Sg{fi) 
if g = 1/G (when the corresponding maximal occupancies coincide). 

First, notice that the function Si/g(/u) cannot coincide with sdp) identically for generic 
G since Sg{fi) does not possess a symmetry like ( pO| ) if 51 7^ 0, 1. Next, we can compare the 



mid-point value Si/a{G/2) = ^^ii ln(G -|- 1) — ^ InG with the corresponding value of ScifJ') 

(26) 



given by (19). Taking into account the inequalities 

tint > (t - l)ln{t + 1) for t>l 



tint < (t - l)ln(t + 1) for 0<t<l, 

we establish that Sg{G/2) > Si/g(G/2) for G > 1. Actually, numerical computations for 
various values of G>1 show that ^^(/u) majorizes Si/G(/f^) everywhere at 0<fx<G, 

SGi^J') > Sr/aifJ-) for G>1. (27) 

For illustration, the case of G = 2 is presented in Fig. 1. 

Thus, we see that for generic G the difference between s^/dtA and Sg(/x) is not small if 
/X is not close to ;U = or /x = G. It is instructive to compare these functions also near the 
end points (in particular, this will provide an additional support to the assertion (^7|)). 

In a vicinity of ;U = equation ( p^ ) is solved as x = ^ + {25g^i — 1)/^^ + 0{p?), and we 
find that Sg(^) = (1 - ln/i)/x -|- (1/2 - (^g,i)/^^ + 0{iJ?), where 

^m,n stands for the Kronecker 

symbol. Comparing this expansion with Sg{n) = (1 — + (1/2 — g)y? + 0{^^) which 

follows from (|^), we infer that for small fi functions Sg{fi) and Sg(/") take close values (up 
to the order jj?) for any choice of G and g. Actually, employing ([T0|), it is easy to derive 
the following expansion for the entropy density of an arbitrary extensive statistics with f{t) 
satisfying (||): 

s{fi) = {Pi- In i,)fi + 0{i^^). (28) 

Thus, for small /i in the thermodynamic limit, the Haldane-Wu and Gentile statistics are close 
not only to each other but to any extensive statistics for which two first Taylor coefficients 
of /(t) are Pq = Pi = 1. For instance, for small /i the Haldane-Wu and Gentile statistics are 



close also to the statistics of 7-ons (18) 



Consider now fi = G — e for small positive e. Then, employing ( |2q ) and the symmetry 
(2C), we derive: Sg{G — e) = (1 — In e) e + O(e^). On the other hand, setting fi = 1/g — e in 



(22), we obtain Sg{l/g — e) = (1 — 2 In (7 — In e) e + O(e^). We see that in a vicinity of fi = G 
the difference between Sg(/x) and Si/aifJ') is not even of first order in e but of order elne (so 
that the corresponding plots are visibly different here, see Fig. 1). 
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Summarizing, it appears plausible that the Gentile statistics of order G > 1 majorizes 
the Haldane-Wu statistics with parameter g = 1/G and these statistics are not close except 
for small values of ii = n/N . 



4 Thermodynamic Bethe ansatz equation 

Now we will consider a relativistic multi-particle system of / different species of particles 
confined to a finite region of the size L. We denote by the number of particles and by Na 
the dimension of the Fock-space related to the species a. In the thermodynamic limit both 
the dimensions of the Fock-space and the size L of the region were the system is confined 



approach infinity but the ratios Ua/L remain finite |12]. For the fraction of particles with 



rapidities between 6 — 1^9/2 and it is convenient to introduce densities 

^Na = Pa{e)ML, Ana = p'a{0) A6 L . (29) 

As usual, the rapidity 9 parameterizes the two-momentum P = m (cosh0, sinh^) and the 
total energy of the system is given by 

' POO 

E[p''] = lY, dO pI{9) ma cosh 9 . (30) 
According to (^) the variable /x in (12) is now related to the particle densities 

Pa{e)=pl{9)/pa{9). (31) 

Thus, we can express the entropy as the following functional 

^ rOQ ^ POO 

S[p,p']=kLY^ / d9 paie)sipa) = klY, / d9{pai9)lnfaiXam-P:m^Xa{0)}, 

(32) 

where each Xa{9) satisfies (p^) and hence it is a function of p'^^{9)/ pa{9). The subscript of fa 
means that different species may have different single-state partition functions. 

In the formulation of the TBA equation for ordinary statistics it is common to introduce 
the so-called pseudo -energies €a{9) such that bosonic (upper sign) and fermionic (lower sign) 
distributions are Pa{9) / Pa{9) = (exp(ea(6')) =F 1)^ . In view of (|lO|), it is natural to define 
the pseudo-energies for an arbitrary extensive statistics as follows 

eai9) = - In xa{9), (33) 

so that p'^a{9) / pa{9) = —9^^(51) In /a(e~'^"^^^). For instance, the distribution for the Gentile 
statistics acquires the following form in terms of the pseudo-energies (with a slightly different 
meaning it was introduced by Gentile for an ideal gas) 

Pm ^ TM=ik^^v{{Ga-k)ea{9)} 

P-^^) Y!i=o^Mkea{e)} 

According to the fundamental principles of thermodynamics the equilibrium state of a 
system is found by minimizing the free energy F. Hence, keeping the temperature constant 
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we obtain the equilibrium condition by minimizing F [p, p*" ] = E [p^ ] — TS [p, p'' ] with respect 
to p^. The equihbrium condition reads 



5pl 5pl bpl t^^^Pb^Pl ■ ^ ' 

The admissible momenta in the system are restricted by the boundary conditions. To 
derive the corresponding quantization equations one takes a particle in the multi-particle 



wave function on a trip through the whole system [12, 13]. The particle will scatter with all 



other particles, which yields the following set of equations determining admissible rapidities 

N 

ex.p{iLmi sinh6'j) Y[ Sij{9i — Oj) = Ki , i = 1, . . . , N , (36) 

where N = J2a^a, and Sij{6) is a two-particle scattering matrix. The constant phases Kj on 
the r.h.s. of ( p6| ) may depend on the particle's species; for our purposes their exact values 
are irrelevant. Taking the logarithmic derivative of (|3^ ) and employing densities as in (ps|), 
we obtain 

I 

ma cosh + 27r ^ {ipab * pim = 27r p,(0) . (37) 
b=i 

Here {u * v){6) = l/(27r) / d6'u{9 — 6')v{9') stands for the convolution and we introduced as 
usual the notation ^pab{0) = — ^^^(In 5af)(^)). 

Employing equations (|To|), (|l2|) and (|3l|), we compute variations of the entropy (^): 

5S/5pl = -\nxa{e) , 6S/6pa = lnfa{xa{9)) . (38) 



Substitution of these relations into the equilibrium condition ( |3q ) together with ( p7D yields 
the desired thermodynamic Bethe ansatz equations for a system in which statistical interac- 
tion between particles of a-th species is described by a single-state partition function fa, 

1 ' 

—ma cosh e = -lnXa{e) + J2 i^ab * In fbixb))i9) . (39) 

'^'^ 6=1 

In general, this set of integral equations cannot be solved analytically. However, one 
can try to solve these TEA equations numerically with the help of the iteration method (as 
we will do for some examples in Section 7). For the ordinary statistics and the Haldane- 
Wu case it is known that such an iterative procedure converges provided that fabi(^) falls 
off sufficiently fast. The most complete proof of this assertion is presented in p^ , p5[ and 
is based on the fixed point theorem. Presumably this proof extends to the case of (^) 
(although some further restrictions on fa{t) may be required), but we will not discuss this 
here. Let us notice only that if (^) can be solved iteratively and <Pafe(^) is symmetric in 6 
(which is always the case for a unitary S-matrix), then Xa{0) is also symmetric in 9. 

We now substitute the expressions for the total energy (^) and the entropy ( |3^ ) together 



with equations (37) and (|3^) back into the expression for the free energy, and obtain 

T ^ POO 

F{T) = -—- J2 ^ria / d0 coshe ln/,(x„(0)) . (40) 



a=l 
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The relation between the free energy and the finite-size scaling function is well-known to be 
c(r) = —QF{T)/{-kLT'^) ||2^. As usual we introduce a variable r = 1/T and set now the 
Boltzmann constant to be one. Then the finite-size scaling function acquires the form 



c(r) 



a=l 



d^ cosh^ ln/,(xa(^)), 



(41) 



Here we assumed that Xa{9) =Xa{—0). Once more for f{t) = 1 + t and f{t) = 1/(1 — t) we 
recover the well-known expressions for the fermionic and bosonic scaling functions. 



It is instructive to compare the above TEA equation (39) and the finite-size scaling 
function (41) with those for the Haldane-Wu statistics |23]: 



1 ' 
—ma cosh e = ln(l + ya{e)) + ^ {^ab * HI + y^')m , 

'^'^ 6=1 



(42) 



c„[r 



E 



d9 coshe ln{l + y-\9)) 



(43) 



Here ^ab{9) = (Pab{9) — '^'^9ab^{9), and g is a matrix which appears on the r.h.s. of (22) if we 
consider a system of several species. As was discussed in [23|, the TEA equation (42) leads to 
an equivalence principle: two systems having the same mass spectra and identical quantities 
^ab{9) are thermodynamically equivalent (as seen from (^), their finite-size scaling functions 
coincide). Our TEA equation (^) does not possess such a feature. More precisely, for (^) 
there exists no way to compensate a difference in statistics by a change of fabid) independent 
of the explicit form of the S-matrix. 



5 Y-systems 

For some class of models it is possible to carry out certain manipulations on the TEA 
equations such that the original integral TEA equations acquire the form of a set of functional 



equations in new variables Ya |14]. These functional equations (commonly referred to as 
Y-systems) have the further virtue that unlike the original TEA equations they do not involve 
the mass spectrum. In the y-variables certain periodicities in the rapidities are exhibited 
more clearly. These periodicities may then be utilized in order to express the quantity Ya{9) 
as a Fourier series, which in turn is useful to find solution of the TEA equations and expand 
the scaling function as a power series in the scaling parameter r. We will now demonstrate 
that similar equations may be derived for a multi-particle system in which the dynamical 
scattering is governed by the scattering matrix related to the ADE-affine Toda field theories 
and the statistical interaction is of a general type considered above. 

Consider the minimal part (i.e., independent of the coupling constant) of the scattering 
matrix of the ADE-affine Toda field theories. As was shown in |27|, upon appropriate choice 
of CDD-ambiguities, these S-matrices satisfy the identity 

I 

e + -) + ^ab[9 - -) = Y^ lac^cbiO) - 2TrIab5{9) , (44) 

c=l 
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where h denotes the Coxeter number, I the rank and / the incidence matrix of the corre- 
sponding Lie algebra. It is then straightforward to derive the "Y-system" 

+ f )ya(^-f) =n(z.(nW))'", (45) 
b=l 

where Ya{9) = x^^{9) and 

Za{t) = t fa{t-^) . (46) 

Equations (|45| ) follow upon first adding (|39| ) at 9 + ^ and 9 — ^ and subtracting / times 



(3£) at 9 from the sum. Thereafter we employ the fact that the masses of an ADE 
affine Toda field theory are proportional to the Perron- Frobenius vector of the corresponding 
Cartan matrix, i.e., Yl,b^ahi^h = 4sin^(7r/(2/i))ma. Then, with the help of (|4^), equations 
dH) follow. 

In comparison with (^) equations ( p5| ) have already the virtue that they are simple 
functional equations and do not involve the mass spectrum. For the Gentile statistics we 
have z{t) = t + l + t~^ + . . . H-t^"^. In particular, for G = 1 we have z{Y) = 1+Y and recover 
the known fermionic Y-system. The bosonic case (G = oo) corresponds to z{Y) = Y'^ / {Y —1) . 
And the case of G = 2 seems to be particularly interesting since here z{Y) = 1 + y + 
possesses an additional symmetry z{Y) = ziY^"^) or, equivalently, z{ea{9)) = z{—ea{9)). 

6 Ultraviolet limit 

In the small r limit the scaling function c(r) becomes the effective central charge of a confor- 
mal field theory [^6| which describes the ultraviolet limit of a given massive model. That is, 
lim,._^o c{r) = Cefj = c— 24/i', where c is the conformal anomaly and h' is the lowest conformal 
weight in the corresponding conformal field theory. 

In order to evaluate c(r) in the r — > limit we substitute the derivative of (39) into (^), 
replacing in both equations cosh 9 by \e.^ ■ Then we integrate the resulting equation by parts, 
assuming that (pab{9) is symmetric and falls off sufficiently fast (typically, (pab{9) = 0{e^^) 
as — > oo). Next we substitute back equation (^). Finally we change variables from 9 to 
x{9), taking into account that 2;(oo) = (as follows from (p9|)), and again integrate by parts. 
All this yields 

Ceff = ^ Ca , (47) 
a=l 

1 , , „ . . r""" dt 

2 

Here we denoted Xa = Xa{9 = 0). These quantities can be determined from the constant 
TBA equations that follow from ( p9| ) when r — > (this derivation requires that Xa{9) become 
constant in a region of 9 of order — In r when r is small; see more detailed explanations in 

/ 

Inxa = -^Nablnfhixb) , a = l,...,/, (49) 
6=1 

where 27r Nab = - d9ipab{9). 

As a simple example we consider the statistics of 7-ons ( [l8[) , assuming for simplicity that 

7 is the same for all species (generalization to different 7a is obvious). Here the integral in 



lnxalnfa{xa)+ -ln/a(t). (48) 
JO t 
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di^ ) can be computed explicitly in terms of the Euler dilogarithm which is defined for t < 1 
as Li2{t) = J2k>it''/^^- We have 

t -y Jo t 7 ^1 «^ 7 

and hence ^ ^ 

Ca = Li2{-^Xa) - — lnj;aln(l + 'yxa) , (50) 

7 ^7 

where Xa satisfies (^9|). Since Li2{-t) = ^Li2{t'^) - Li2{t), we can rewrite ( ^OD for positive 
7 employing the Rogers dilogarithm which is defined as L{t) = Li2{t) + ^lntln(l — t) for 
0<t<l. Then, using the Abel identity L{f) = 2L{t) - 2L{j^^), we obtain 

Ca = -L(-^^) + ^lnjHl + ^Xa). (51) 
7 Vl + 7x„/ 27 

Notice that this expression does not diverge for small 7; lim^_^oCa = Xa{l — ^Inxa) is the 
value corresponding to the Boltzmann statistics. 

Consider now the Gentile statistics, assuming that the order G is the same for all species. 
In this case we also can compute the integral in ( ^8[) explicitly: 

X 



-In fait) = I -In = L 12 0^ " T?" T = ^^^(x) 



and we see that for the Gentile statistics expression ( p8| ) acquires the following nice form 
involving only the Rogers dilogarithms: 

c, = L(xJ-^l(x^+i) . (52) 



Introducing, in agreement with (|33D, the constants = — Inx^, we obtain from (52) for 
G = 00 and G = 1 

Ca = L{e''^) and = ^(6-'") - iL(e"2'") . (53) 

These are the well-known values for the bosonic and fermionic statistics Q. Notice that in 
the latter case one usually derives Ca = L{j^^). This is equivalent to our formula ( |5^ ) due 
to the Abel identity. In fact, the Abel identity can be used in a similar way for any G of the 
form G = 2*" — 1 . In this case we can rewrite ( |52|) as follows 

" 1 ^ (r . 



k=l 



In this form Ca is a manifestly positive and monotonic function of Xa (the dilogarithm L(t) 
grows monotonically for < t < 1). Actually, this property holds not only for the Gentile 
statistics but in general case as well. To prove this statement we take a derivative of (^) 
and use (0) and the formula (^) for the entropy density. This yields 

a^^^^^dJ^. (55) 
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As we discussed in Section 2, we should consider only such models where s{fj,{x)) >0. In this 
case the r.h.s. of ( |55[ ) is positive. Hence, Ca is always a monotonically increasing function of 
Xa (and therefore a monotonically decreasing function of ea). Moreover, taking into account 
we infer from (^) that Ca(0) = 0, thus we also proved positivity of Ca{xa)- Since these 
properties are crucial for a physical interpretation of Ca, let us underline that they hold only 
if the first condition in (0) is fulfilled. Let us notice also that the properties of Cq together 
with ( ^9| ) imply that for any model such that Nab ^ (which is the case, for instance, for all 
the ADE-affine Toda models Q) the value of Ca does not exceed that of the corresponding 
free model (where all Nab = and consequently all Xa = !)■ Hence, ( |4^ evaluated for Xa = 1 
yields the upper bound on Ca, 

ca< / -In fait). (56) 
JO t 

Notice, however, that for a statistics of the type HIb with i? < 1 we have Xa < 1, so that 
in this case the upper bound on Ca will be lower than ( j56| ) . From a physical point of view 
this implies that such a statistics cannot emerge in models with weak interaction and, in 
particular, in a free model. 

It was argued in Q that the quantity :^Ca (for the ordinary statistics) can be interpreted 
as a massless degree of freedom associated to the a-th species in the massive model. The 
properties of Ca which we proved above show that the interpretation of :^Ca can be extended 
to the case of a generalized extensive statistics as well. In this context, it is important to 
remark that for the ordinary statistics (actually, for the Gentile statistics of any order) we 
have ;;^Ca < 1 as follows from ( |52[ ) (recall that L(l) = Li2(l) = tt^/G) if all Nab are non- 
negative. However, for a general statistics the upper bound ( |56D on :^Ca depends on the 
choice of fa and may not equal to one. For instance, taking fa{t) = (1 — t)~™' with m > 1, 
we obtain -^Ca = m for the corresponding free theory. In such a case, the massless degree 
of freedom of the corresponding free particle exceeds that of a free boson. This can possibly 

2 

imply that we have to restrict the choice of statistics to such fa that Ca < 



7 Examples related to affine Toda models 

7.1 Ising and Klein- Gordon models 

The most simple examples which illustrate the features outlined above more concretely are 
the Ising model (Ai-minimal affine Toda field theory) with S{6) = — 1 and the Klein-Gordon 
model with S{6) = 1. In both cases equation (^) is solved trivially, x{9) = exp{— rm cosh 0}. 
Then for the Gentile statistics (|^) we can compute the entire scaling function (p]): 

Ca{r) = ^^°°decosh0(ln(l-e-(^+^)™"°'^^)-ln(l-e-^"^^°'^^)) 

= ^f]y(Kiikrm)-Ki{{G + l)krm)) , (57) 

where Ki{t) is the modified Bessel function. We depict Cc{r) in Fig. 2, referring to them by 
a slight abuse of notation as Ai. One observes that the difference in statistics affects most 
severely the ultraviolet region. The scaling functions for different G converge relatively fast 
towards each other in the infrared regime. These appear to be common features of a general 
deformed statistics. 
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As seen from Fig. 2, the behaviour of C(j(r) in the ultraviolet region (there is no first order 
term in the small r expansion of Cc{r)) is similar for any finite G to that of the fermionic 
scaling function and is different from the bosonic case (where the first order term is present). 
To explain this we observe that, as follows from ([57|), 

CG(r)=Cb(r)-^Cb((G + l)r), (58) 

where c\,{r) is the bosonic scaling function. Now it is obvious that the first order terms on 
the r.h.s. of (|58|) cancel each other for any finite G. 

The well-known property of the asymptotic behaviour of the modified Bessel function, 
limt^otKi{t) = 1, allows us to derive from ( |57| ) 

fc=i 

which is in agreement with ( |5^ since Xa = 1 now. It is worth noticing that c'^ can be 
identified for every positive integer G as the effective central charge of a minimal conformal 
model M.{s,t) (such that st = 6{G + 1)). In particular, G = 1,2,3,4 correspond to the 
M{3,A), M{2,9), M{3,8), M{5,6) minimal models, respectively. 

For the 7-ons ( p!8| ) the scaling function related to the Ai-models can also be computed: 

cJr) = ^ de cos\ie\n{l+-ie-''^''°''^% = —^y -{--lf'^Ki{krm). (60) 

Using again the asymptotic behaviour of the modified Bessel function, we obtain c'^ = 
\\m.r^QC^{r) = —^^Li2{—^) in agreement with (^). Here, unlike in the Gentile case, it is 
not easy to see whether the effective central charge has rational value for a given 7. As 
we discussed in Section 2, positive 7 is to be of the form 7 = 1/d with d positive integer. It 
appears then that the only d leading to rational value of is d = 1. On the other hand, 7 
varies continuously in the range — 1 < 7 < 0, so that any rational value of not exceeding 
^ occurs here. For the Boltzmann statistics Cg^ = \s irrational. Finally, 7< — 1 is the 
case of a type Illb statistics with i? < 1, which, as we discussed above, cannot emerge in a 
free model. 



7.2 Scaling Potts and Yang-Lee models 

Next we consider the scaling Potts model, which was studied previously in the TBA frame- 
work for the fermionic statistics in |[l^ and for the Haldane-Wu statistics in |23]. The two 



particles in the model are conjugate to each other and consequently their masses are the 
same, mi = 1712 = m. The conjugate particle occurs as a bound state when two particles 
of the same species scatter. The S-matrix of the scaling Potts model equals the minimal 
S-matrix of the ^2-affine Toda field theory. The corresponding quantities (Pabi^) that enter 
the TBA equation are given by 

Consider the scaling Potts model with both species of particles obeying the Gentile 
statistics of the same order G. In this case the Z2-symmetry of the model is preserved so 
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that xi{9) = X2{0) = x{6) and (|39|) reduces to a single integral equation. It appears to be 
impossible to find an analytic solution x{6) to this equation, but it is straightforward to solve 
it numerically. Taking x^^\0) = exp{— rm cosh 0} as the first approximation, we can iterate 
( p9| ) as follows 

Convergence is achieved relatively quickly (it depends on the value of mr). The results for 
G = 1 (which is in complete agreement with the calculation in [jl^ ) and for G = 2 are shown in 
Fig. 3. To make contact with the literature, we introduced the quantity L{9) = In fG{x{6)). 
One observes the typical plateau L{9) = const for some region of 6 when mr — > 0, which 
is required to derive (^). We see also that for the same value of mr the functions L{6) 
corresponding to different G have similar profiles but different heights of the plateau. 

Having solved the generalized TBA equation, we can substitute x{6) into (41) and com- 
pute the entire scaling function cdmr). The result of the numerical computation is shown 
in Fig. 4. The behaviour of the scaling functions is similar to what we have already observed 
in the Ai case. 

Due to the fact that the scattering matrices of the scaling Potts model and the scaling 
Yang-Lee model are related as S^^{0) = S^^{9)S^2{^)^ the generalized TBA equations for 
the two models coincide. The only difference is that the scaling Yang-Lee model has only 
one particle. Therefore its scaling function equals a half of that for the scaling Potts model. 

7.3 Comparison with Haldane-Wu statistics 

In Section 3 we have seen that the entropy density corresponding to the Gentile statistics 
majorizes that for the Haldane-Wu statistics. It appears (at least for certain models) that 
an analogous relation holds for the corresponding finite-size scaling functions. Specifically, 
if Ccir) is the scaling function ( pT|) of Ai or A2 minimal affine Toda model where all species 
obey order G Gentile statistics and Cg{r) is the scaling function ( ^3[ ) of the same model where 
the Haldane-Wu statistical interaction is of the form Qah = ^^ab, then for all r > we have 

Coir) > Cg{r) , where gab = ^Sab and G > 1 . (63) 

For the Ai case we can rigorously prove relation (|63| ) since the solutions of the generalized 
TBA equations (EdI) and (E2|) are found explicitly: 



g lny{e) + {I - g) ln(l + y{e)) = mr cosh 6 = - lnx{0) . (64) 

This equation allows us to express x{9) in terms of y{0): 

x{e)=y'^9){l + y{e)r-\ (65) 

Moreover, since mr cosh 0>O, we infer from ( |64|) that 

y{e)>yo, where y9 = {1 + yof'^ (66) 

As seen from (^) and (|^), in order to establish (63) it suffices to show that fc^i/g{x{9)) > 
1 + y~^{9). The latter relation is equivalent to the inequality 

1 - (1 + yie))-^ > y3{9) (1 + y{e))-^ , (67) 
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as can be verified by simple algebraic manipulations using formula (^). Now it is elementary 
to check that ( |67| ) is valid provided that (^6|) holds. Thus, we proved the assertion (|^) in 
the Ai case. For illustration, the scaling function Cg{r) for g = l/2 Haldane-Wu statistics is 
shown in Fig. 2. It lies everywhere below Ca{r) for G = 2 Gentile statistics with the maximal 
difference reached in the ultraviolet region. 

For a higher rank case, already for A2, the generalized TEA equations are integral and 
it is not clear whether we can compare Xa{9) and ya{G) analytically. However, numerical 
computations for various G show that (^) holds in the A2 case as well. For illustration, we 
present in Fig. 4 the scaling function Cg{r) for the Haldane-Wu statistics with gab = \^ab- 

A particular consequence of (|63| ) is the following inequality for the corresponding central 
charges: 

cf>cf, for gab = h^ab and G>1. (68) 

Here d^^ is found from (47), (49) and ( ^2| ) whereas is the effective central charge related 
to the Haldane-Wu statistics with gab = -^^ab and is given by [23| 



= 11^(7^:7) ' ^^6^6 Hl+ya)=Y,{Nab + 95ab)H'^ + ya^)- (69) 



eff ^ 6 



Again, in the Ai case we can prove relation (^) directly. Here c^^ has a simple form (^) 
and therefore (|^) reduces to 

' >4l(^1, (70) 



1 + 5- vr^ ^l + yo' 

where yo is defined as in (|6^). Notice that g < y^ < 1 for (7 7^ 0, 1 (since yo < g would 
contradict ([2^) ). Therefore, ^ < < Using now a property of the dilogarithm, 

2 2 
L{t) < y t for i < t < 1 ; L(t) > for < t < i , (71) 

we establish the desired inequality (|70|). 

As we discussed above, the order G of the Gentile statistics is to be a positive integer. 
However, if we regard (^7|), (|4^ ) and (|52|) as a definition of the quantity c^''^, then (70) extends 
to any G > 1 (since the proof does not require that G be an integer) as shown in Fig. 5. 
Furthermore, considering G in the range 0<G< 1, we can prove with the help of (p6[) that 



here yo < g and hence :^L{j^^) > due to (|7l|). Thus, for this range of G we have a 
reverse inequality, i.e., 



<cf, for gab = h^ab and 0<G<1. (72) 

As we will see below, this inequality holds for higher rank cases as well. But let us stress 
agam that cf in (|72|) is a formally defined quantity which is not related directly to the 
Gentile statistics. 

The inequality for the central charges ( |68D is a weaker statement than the inequality for 
the scaling functions (|6^). Furthermore, (|6^ ) involves dilogarithms, which makes its proof 
more complicated (already in the Ai case we used a rather subtle property (|7l|) ). But the 
advantage of (|68|) for a numerical verification is that it relates numbers and not functions 



c, 



as (63). Moreover, since c(r) is a continuous function, validity of (pq ) implies that 



holds at least in some ultraviolet region. Therefore, we will discuss below validity of (|68|) for 
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several affine Toda models. As a mathematical by-product, this provides us with interesting 
dilogarithm inequalities. 

In the A2 case, eq. (^) contains dilogarithms in a more involved way. As seen from (^), 



(52) and (p9[), it is equivalent to the following inequality 

L(xo) - -^l(xI~^^) > , for < 5 < 1 . (73) 

Here xq and yo are determined from the following equations that follow from ( ^9| ) and (|69|) 
upon noticing that A^ii+A''i2 = l 

xofi/a{xo) = l and y^+f = (1 + yo)^ . (74) 

It is worth remarking that, as was noticed in ||2^, ^L( ^^-^ ) coincides with the effective 
central charge of the Calogero-Sutherland model with the coupling constant X = g. 

Numerical computations show that ( [73|) indeed holds for any < g <1, that is for any 
G>1 (here again we need not to restrict G to be an integer). Furthermore, like in the Ai 
case, inequality ( |73| ) reverses for < G < 1, i.e., for g > 1. For illustration, we depict the 
corresponding and Cg^ in Fig. 5. This result provides an additional support to our claim 
that in the A2 case the entire scaling function Ca{r) majorizes Cg{r) for G>1. 

It is interesting to remark that, as follows from ( ^2|) and (|69|), the bosonic version (G = oo 
or, equivalently, g = 0) of the j42-minimal affine Toda model has Ca = |. That is, then each 
species in the scaling Potts model (or the only species in the scaling Yang-Lee model) appears 
to be a free fermion, whereas the entire central charge is that of a free boson (see also related 
comments in p^]). That is why all curves on Fig. 5 converge to the same value c=l. 



7.4 Higher rank cases 

Now we discuss briefly the simplest minimal affine Toda models which have particle species 
with different masses. We will compare numerical values of and dg^ again setting gab = 
^Sab- For this purpose we use ( ^9| ) and (|69[), taking into account that the A^-matrix in 
these equations is given by A'g = /g(2 — Ig)^^, where /g stands for the incidence matrix of 
the corresponding Lie algebra g. 

In the ^3 case the masses of the three species are mi = 777,3 = m2/\/2- Then, computing 



cf by formulae (49) and (52), we find: c^** 0.77, c^** = 1, c*^" « 1.12, c<=" « 1.16 for 



G = |,l,2,cxD, respectively. The values of Cg^ corresponding to the Haldane-Wu statistics 
with gab = ^Sab were found in @]: d"^ ^ 0.89, c"^ = 1, d'^ ^ 1.07, d"^ « 1.16 for 



g = 2,1, i,0, respectively. Carring out analogous computations for the A4 case, where the 
masses of the four species are 7771=7774, 7772 = 70,3 = 7771 (\/5-|-l)/2, we find (Gentile statistics): 
d"^ ^ 0.92, d"^ = 8/7, c^^ « 1.25, d'^ « 1.28 for G = 1, 1,2, 00, respectively; and (Haldane- 
Wu statistics): c"^ 1.05, c^*^ = 8/7, c"^ ^ 1.20, c^^ « 1.28 for g = 2,l, i,0, respectively. 

(2) 

One can also consider a folded algebra case with two different masses, namely A\ , where 
77i2 = 777i(\/5+l)/2. In general, one assigns to the A^2rl minimal affine Toda model a tad-pole 
graph ^ such that the corresponding incidence matrix differs from that of An only by 
1 in the lower-right entry. Then one can check that Xa in ( p9D and Ua in (^) coincide with 
those for the A2n case (a well-known fact for the fermionic statistics O]). Consequently, 
for any statistics we have c (A2n ) ~ 2*^ (^2n)- Therefore, the data for A\ follow from the 
above results for A4. 
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Finally, in the D4 case, where the masses of the four species are mi = 1713 = 1714 = 1712/ V^, 
we obtain for the Gentile statistics: c""^ ^ 0.82, d"^ = 1, c""^ « 1.08, c"^ « 1.09 for G = 
^,1,2,00, respectively; and for the Haldane-Wu statistics: c^^ « 0.93, d^^ = 1, d^^ « 1.04, 
d^ PS 1.09 for g = 2,1, \, 0, respectively. 

We see that all the above results are in agreement with (^) and (|7^). This, together 
with the more detailed results we obtained in the Ai and A2 cases, allows us to conjecture 
that these inequalities and the inequality (^3|) for the scaling functions hold actually for all 
simply-laced minimal affine Toda models. 

We remark also that in all the computations of d^^ in this subsection we always have the 
rule that the value of Ca is smaller for the species which has heavier particle (this was known 
for the fermionic statistics Q). This observation provides an additional support to the above 
discussed interpretation of -^Ca as a massless degree of freedom of the corresponding particle. 



8 Conclusion 

Summarizing, we considered possible types of a generalized extensive statistics and estab- 
lished properties of the corresponding entropy densities. In particular, we established nu- 
merically (and gave some supporting analytical arguments) that the entropy density for the 
Gentile statistics of order G > 1 majorizes that for the Haldane-Wu statistics if they cor- 
respond to the same maximal occupancy. Further, we derived the thermodynamic Bethe 
ansatz equation and the finite-size scaling function c(r) for a relativistic multi-particle sys- 
tem obeying a generalized extensive exclusion statistics. We put particular emphasis on the 
ultraviolet limit of such a system. We derived an expression for an effective central charge 
in the general case and showed that for the Gentile statistics it acquires an elegant form in- 
volving dilogarithms. We discussed a physical interpretation of the 'partial' central charges 
Ca and argued that it possibly leads to restrictions on the choice of statistics. Finally, we 
observed (and partially proved) that the majorizing properties of the Gentile statistics with 
respect to the Haldane-Wu statistics extend also to finite-size scaling functions and central 
charges related to the minimal simply-laced affine Toda models. 

In the presented analysis of thermodynamic properties of a relativistic multi-particle sys- 
tem the single-state partition function f{t) played a key role. For a given model, if we know 
f{t) (at least in the thermodynamic limit) and the quantities ^ab{G)^ then we can compute, 
at least in principle, the entire scaling function which provides complete thermodynamic 
description of a system. In principle, f{t) should be determinable from the correspond- 
ing Hamiltonian. However, this might be a non-trivial task since the exclusion statistics 
of particles can differ from the exchange statistics of the fields present in the Hamiltonian. 
Therefore, it may be more practical for a given system to fix an exotic statistics a-priory and 
investigate whether it leads to realistic physical properties. In particular, one can obtain 
information about the ultraviolet limit and see if it corresponds to an appropriate confor- 
mal field theory. As we have seen above, the different types (|8|) of f{t) lead to significantly 
different thermodynamic properties. The question, which of these statistics can emerge in 
physical models, remains open. 

Another open problem is to prove rigorously the assertion ( p7| ) that the entropy density 
for the Gentile statistics majorizes that for the Haldane-Wu statistics. Such a proof will 
probably provide a deeper insight into general properties of exotic exclusion statistics. Also 
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it would be interesting to verify the majorizing inequalities for the scaling functions (p^) 
and central charges (^) related to the affine Toda models and to understand whether such 
a majorization is a model independent feature of the two involved statistics. 
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Fig. 1: Entropy densities for G = 2 Gentile statistics (the upper curve), for g = 
statistics (the middle curve) and for the fermionic statistics (the lower curve). 
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Fig. 2: Scaling functions of the Ai-minimal afhne Toda model for the Gentile statistics of order 
G = 1,2,3,4 (solid curves, from down to up), for g = 1/2 Haldane-Wu statistics (the dashed curve) 
and for the bosonic statistics (the dotted curve). 
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Fig. 3: Solution L{6) of the TBA equation for the A2-minimal affine Toda model for the scaling 
parameter mr = 0.1 (dotted curves), mr = 0.01 (dashed curves), and mr = 0.001 (solid curves). For 
the same value of mr the lower curve corresponds to the fermionic statistics and the upper curve to 
G = 2 Gentile statistics. 
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Fig. 4: Scaling functions of the A2-mininial afSne Toda model for the Gentile statistics of order 
G = 1, 2, 3, 4 (solid curves, from down to up) and for the Haldane-Wu statistics with gat = ^Sab (the 
dashed curve). 




G 



Fig. 5: Effective central charge of the Ai and j42-minimal affine Toda models for the Gentile statistics 
of order G (lower and upper solid curves) and for the Haldane-Wu statistics with gab = ^Sab (lower 
and upper dashed curves). 
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